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In models in which neutrinos are light, due to a low scale of symmetry breaking, 
additional light bosons are generically present. We show that the interaction between 

m 

' diffuse supernova relic neutrinos (SRN) and the cosmic background neutrinos, via 

CN ■ exchange of these light scalars, can result in a dramatic change of the supernova 

<N , 

, (SN) neutrinos flux. Measurement of this effect with current or future experiments 

in 
o 



X 



can provide a spectacular direct evidence for the low scale models. We demonstrate 



' how the observation of neutrinos from SN1987A constrains the symmetry breaking 

scale of the above models. We also discuss how current and future experiments may 
confirm or further constrain the above models, either by detecting the "accumulative 



resonance" that diffuse SRN go through or via a large suppression of the flux of 
neutrinos from nearby ^ O (Mpc) SN bursts. 

I. INTRODUCTION 

The observation of neutrino flavor changing from solar fl], [2|, [|, atmospheric Q and 
terrestrial || || neutrino data has provided firm evidence for neutrino flavor conversion. The 
recent new Super-Kamiokande (SK) data on the L/E'-dependence of atmospheric neutrino 
events J/J, L being the distance traveled by neutrinos of energy E, and the new spectrum data 
from terrestrial experiments || [|, has yielded for the first time evidence of the expected 
oscillatory behavior. This strongly favors non-vanishing sub-eV neutrino masses. These 
outstanding developments on the experimental side of neutrino physics have placed a distinct 
burden on theorists — to understand what is the origin of these tiny neutrino masses. 



The most elegant and popular solution to this puzzle is the seesaw mechanism [1C]. In 
this scenario one assumes that lepton number is violated at some high scale Al in the form 
of right-handed neutrino Majorana masses. This induces, at a lower scale, an effective 
operator of the form 0(1) x {LH) 2 / Al , where L denotes a lepton doublet and H the Higgs 
field. The oscillation data then imply that Al ~ 10 14 GeV . While the seesaw mechanism is 
very appealing from the theoretical side, it is unlikely to be subject to direct experimental 
test sometime in the near future. An additional virtue of the seesaw mechanism is that 
it can naturally provide a platform for generating the observed baryon asymmetry of the 



universe through leptogenesis |TT[. Introduction of such a high scale, however, requires a 
mechanism for electroweak-symmetry-breaking-scale stabilization which typically leads to 
various moduli/gravitino problems in the context of cosmology. Thus it is important to 
explore alternate origins for neutrino masses. 

One such alternative is the late neutrino mass framework that induces small neutrino 



masses due to a low scale of symmetry breaking fL2| . This idea points to a completely 
different understanding for the origin of neutrino masses. The neutrino masses are protected 
by some flavor symmetry different from the one related to the charged fermion masses. 
When this symmetry is (say spontaneously) broken by a set of flavor symmetry breaking 
vevs, /, of fields 0, the neutrinos acquire masses from LNH for Dirac neutrinos, or 

(jjf^) LNH + MrNN for Majorana neutrinos, where N denotes a right handed neutrino, 
L, H stand for the SM lepton doublet and Higgs fields respectively and Mp is a scale in 



which flavor dynamics takes place |L3| . We want to stress that these textures do not depend 
on the details of the symmetry mechanisms, whether global |l2j or gauge [ 14| . Furthermore 
a similar scenario can be realized via strongly coupled dynamics where the compositeness 
scale is given by / [[15], [U| . 

With this alternate scenario, it is then of immediate import to delimit the allowed range 
for the symmetry breaking scale, / at which new physics (NP) appears. Since the principal 
consequences of the symmetry breaking are neutrino masses and the relevant new degrees 
of freedom couple only to neutrinos, direct experimental limits on the parameters of this 
model are unlikely to be attained. In fact, the strongest limits on / come from cosmology 



and astrophysics rather than from laboratory data |12|]. As will be discussed shortly, there 
are generically associated with this mechanism some extra light degrees of freedom. In the 
case that the number of these exceed present bounds from big bang nucleosynthesis (BBN), 
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the requirement that these not be in thermal equilibrium during BBN gives a limit on / of 
approximately [12|, [14|, 16, [L7|, [18 

/>10keV. (1) 

A similar bound is obtained by demanding that SN cooling not be modified in the presence 
of the above additional fields. 

It is remarkable that this framework with a low NP scale, / ^ Aewb where Aewb is the 
electroweak symmetry breaking (EWB) scale, cannot be excluded by direct experimental 
data. In many late neutrino mass models, there are degrees of freedom beyond . These 
additional degrees of freedom can yield indirect signals provided that standard cosmology 
is assumed. In the case in which neutrino masses are protected by global or approximate 
symmetries |12| or the case with strong dynamics (in which chiral symmetries are being 
broken by the condensate) [TH], light pseudo-Goldstone bosons (PGB) field are typically 



present. Similarly, in models with gauge symmetries [H the corresponding gauge boson 
masses are suppressed, relative to /, by an additional gauge coupling g, and therefore play 
a role similar to the one played by the PGBs. This additional light fields interact with the 
plasma through their coupling to the neutrinos. This happens even below the BBN phase 
transition and may leave a trace in the observed cosmic microwave background radiation 
(CMBR) jg. 1 

The focus of this work is to investigate other more direct ways of testing the low scale 
models of neutrino masses. We find that such a possibility of a more direct probing of 
this class of models, at present or in the not-too-distant future, does exist. The desired 
signal would consist of a dramatic modification of the supernova neutrino flux (diffuse or 
burst) through interaction between the these neutrinos and the cosmic background neutrinos 
(CBN). These interactions are mediated by the new scalar particles introduced by the NP. 

In Section y we discuss the dominant processes which modify the incoming SRN fluxes. 
We divide this section into two: in |U A| we describe the resonant process which happens only 



in a narrow range of parameter space, but leads to a spectacular signal through what we 
denote as accumulative resonance; in [II B| we consider the non-resonant processes which, in 
conjunction with data from SN1987A, yield a lower bound on / comparable to the one from 
BBN. Also, at the beginning of Section [II] we summarize bounds on the model parameters 



1 see also ]l9| for related analysis. 
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imposed by BBN. These bounds are sketched in Fig.[8], along with the parameter range for 
which the resonant process occurs. Finally, we conclude in Section [ITT . 



II. MAIN IDEA, FORMALISM AND IMPORTANT PROCESSES 

Our main idea in this work is to show that the presence of the additional light bosons, 
required in the late neutrino masses framework, can introduce a significant interaction 
between the SRN and the CBN. This, in some region of the model parameter space, can lead 
to a measurable modification of the incoming SRN flux. The typical SRN energies are above 
average solar neutrino energies and below the atmospheric ones. Consequently, this flux is 



likely to be observed by SK [£(J and KamLAND [21| in the near future, or by successor 
experiments. Thus, there is a window (although not very wide) in which we can observe the 
presence of both this extra light degrees of freedom and the CBN! In this part we introduce 
the relevant part of the Lagrangian and discuss important processes which yield the signal. 

We first discuss resonant processes and present the phenomenon of accumulative 
resonance which can yield a possible signal. Then we move to discuss non-resonant processes, 
which, through the observation of neutrinos from SN1987A yield a bound comparable with 
the BBN one. Other implications of these processes, related to experiments envisioned for 
the near future, are also discussed. 

Below EWB scale and close to the neutrino flavor symmetry breaking scale the effective 
Lagrangian can be written as 

CP V = £ kin + y v( j>uN + V{(j>) , = C kin + y v( j>uu + V{<f>) , (2) 

where £® ,M stands for the Dirac and Majorana case respectively, £km denotes the kinetic part 
(for the gauge case this contains interaction between and the additional gauge bosons 
v represents an active neutrino, V(<j)) is the scalar potential (for the global case this contains 
interaction between and the additional Goldstone bosons jl4j|), and flavor and spinor 
indices are suppressed for simplicity. The above implies that 

m v = yj . (3) 

As we shall show below our signal is similar in both the Dirac and Majorana cases. For 
simplicity through our discussion below we omit the effect of neutrino mixings (apart from 
the discussion related to the SRN flux). 



In order to establish a reference point, we first pause to summarize the bounds on various 
models imposed by BBN constraints in terms of the Yukawa couplings y v . This will enable 
an evaluation of the feasibility of our program. 

1) The minimal model is of Majorana neutrinos with Abelian symmetry. We assume 
that the symmetry breaking scale, /, is below the BBN temperature of about 1 MeV. Then 
during the BBN epoch we cannot separate the Goldstone and the scalar (higgs) as they are 
a single entity, a complex scalar field. The updated BBN bound on the number of neutrinos 
is N = 3.24 ± 1.2 at 95% [22|. The complex scalar adds 8/7 (neutrino) degrees of freedom, 



so this additional degree of freedom can be accommodated with the BBN bound above. 
However, there are other cosmological bounds, such as SN cooling rate. If the BBN bound 
on relativistic degrees of freedom should decrease by a significant amount, then the yukawa 
y v would be subject to the upper bound obtained in the next paragraph. 

2) In the non-Abelian Majorana models, typically several complex scalars are present, 
which are not permitted to be by BBN considerations. Thus, in this case y u must be 



bounded from above to ensure decoupling. This bound was derived in Ref. |1| and we have 
also calculated this bound (as a check) by considering all the processes that would produce 
G's. Recoupling via the 2 — > 1 process vv — > G takes place as the temperature falls to some 
value T rec determined by equating the decay rate at T rec to the Hubble expansion rate: 



M G y 2 v M G 8n^g 

3T rec 16tt V 45 M 2 Pl ' U 
where g is the number of degrees of freedom at T rec . By requiring T rec < T B bn we find 

y u < 6 x 10- 7 (keV/M G ) (5) 

3) Finally, for the Dirac case, the absence of a negligible population of right-handed 
(sterile) neutrinos (N) in the bath disallows the reaction vN — > G, so that G's can only be 
produced via vl^l — *• G G (via t channel N exchange). Requiring that this process be out 
of equilibrium at Tbbn yields a BBN bound of 

y v < 1 x 1(T 5 . (6) 

The s-channel process requires a chirality flip which makes the bound weaker, as pointed 



out in Ref. [12]. Note that this bound in independent of the Goldstone mass. 
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A. 



Resonance, Accumulative Resonance 



The simplest possible process which will modify the spectrum is the resonant production 
of one of the above light bosons in the collision of an SRN and a CBN. For simplicity 
we shall assume that the boson couples only to neutrinos with a strength y v (this is a 
good approximation for the case in which y u < as discussed below). Thus after being 
produced the boson, say a scalar, will decay back yet to a pair of neutrinos where in our 
frame they have an energy spectrum flatly distributed between and the resonant lab 
energy E*. The diagram describing this process is shown in Fig. [I]. For clarity, we frame 



FIG. 1: Diagram representing resonant scattering. G denotes a pseudo-Goldstone or a gauge boson. 

the discussion which follows in terms of Dirac neutrinos. Except for differing dynamics in 
the high temperature environment of the supernova, the Majorana case is similar. 

The resonant scattering affects the incoming SRN spectrum: the energy of the incoming 
SRN is now divided between the two decay products, so that we expect to observe a depletion 
in the expected spectrum for incoming neutrino with the appropriate energies. In addition, 
since this process is effectively 1— »2, one may expect some depletion in the observed incoming 
neutrino spectrum at appropriate energies. 

Two questions are in order: 

1. In view of the small upper limit on the interaction strength (|^) and the low density of 
the background neutrinos, will the resonance process indeed produce the depletion discussed 
in the preceding paragraph? 

2. If the answer to (1) is affirmative, can the resultant depletion be observed by present 
or near future experiments? 

Both these questions will be answered in the affirmative. 
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1. Resonance: no cosmological expansion 



First, we consider the case of no expansion and begin by estimating the mean free path 
(m.f.p) ARes for the resonance process. 

The cross section, written in Breit-Wigner form for the process in Fig. [I] is roughly given 

by 

aRes ~ 16^ {Ml - sf + MlYl ' (?) 

where G stands for the gauge/ Goldstone bosons and s is the square of the center of mass 
energy. In addition Y v is the decay width of the boson into neutrino pair, 

v 2 Mr 



Consider the case in which the SRN energy, E, is on resonance, E 



* 2 



so that 



7T 



(10) 



Consequently on the resonance 3 the m.f.p is give by 



ARes- V-~4> (11) 

where n v = 37rr(3)C(3)T i 3 ~ 56(1 + z) 3 cm -3 , n° is the present background neutrino density 
and T v is the background neutrino temperature with T° ~ 1.6 x 10~ 4 eV being the present 
one. On carrying out the numerical substitution and using Eq. ^| we find 4 

x 1 2t7vE* 6 m v E* 

ARes ~ ~ 7T~ ~ 5 X 10 PC . (12) 

nOcr Res 7m° u * 5 x 10" 2 eV lOMeV V ; 

Since this is much smaller than a typical distance traveled by a SRN neutrino, we find that 
the answer to Question 1. is positive for practically any value of y v . That is, if a neutrino 
is produced with the appropriate resonance energy then the process will go through. 



2 We shall neglect here the effect of thermal broadening, effectively assuming that the background neutrinos 
are at rest. We further discuss this point below. 

3 The result is qualitatively the same if one averages over the width of the resonance. 

4 Note that for the lighter neutrino the m.f.p is even somewhat shorter. 
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However, this analysis also implies that the answer to the second question is negative. 
The extent of the dip in the SRN flux is controlled by the width of the boson G . This width 
is tiny (Rl), rendering it impossible at present or in the near future to detect such a narrow 
depletion in the SRN flux. The results are markedly different when cosmological expansion 
is included, so that we now turn to consider this case. 



2. Resonance: case of cosmological expansion 



We start by finding the conditions on the coupling for which there is sizable resonant 
degradation of of the original (not products of G-decay) flux of supernova neutrinos. 

The the probability P(E, z) that a neutrino, created at red shift z, with energy (1 + z)E 
arrives unscattered at the detector with energy E is given |23| by an integration over proper 
time (converted to an integral over intermediate red shifts z): 

dz 



P(E, z) = exp 



— n u a uv ^{2m v {\ + z)E) 



(13) 



Jo H(z)(l + z) 

where a is the resonant scattering cross section, a function only of s = 2m u (l + z)E, the 
(cm. energy) 2 at the time of scattering, n u is the neutrino density (per flavor) at redshift 
z and H(z) is the Hubble constant at redshift z. 

For the purposes of this section, it is sufficient to employ a 5-function approximation for 
the cross section, 

7T 



o = -yi S(s-M^) 



(14) 



so that the integral in (|T^) is trivially done, with the result 



P(E, z) = exp 
for 



A M G 



E 



H(E*/E) 



E* 



l + z) 



<E <E* 



(15) 



and P = 1 otherwise. Here E* is the resonance energy Ml/{2m v ). There will be large 
depletion of the initial SRN flux in this entire domain if 



2 < 



AM 2 G H 

which gives (for Hq = 70 km s" 1 Mpc" 1 ) 

Vu £ 4 x 10' 



M, 



G 



1 keV 



(16) 



(17) 
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In the Dirac or Abelian Majorana cases, this permits a reasonable window of more than two 
orders of magnitude in y v (see Eq. ([])) for substantial depletion of the SRN flux. For the 
non-abelian Majorana case, the window narrows to a bit more than an order of magnitude 
(Eq. (|^)). In this "strong coupling" regime we take P = in the domain in Eq. (|15|). Thus, 
with the appropriate constraint on y u , Question 1. is answered in the affirmative: there will 
be substantial depletion of the original flux due to resonant scattering. We now proceed 
to see how the cosmological evolution permits a signal to be formulated for the resonant 
scattering. 



3. Accumulative resonance 

In the domain (|15|), there will be resonant absorption out of the original neutrino flux, 
but some replenishment as well, from neutrinos re-emitted in the decay of a G produced 
in the domain in Eq. (pT5|). More specifically, suppose that a neutrino emitted with energy 
e > E* from a source at redshift z undergoes resonant scattering at redshift z < z, so that 

£* = e— . (18) 
l + z 

This is followed by the emission of a decay neutrino with energy E' = fE* , < / < 1 
immediately following emission. The flatness of the emitted-neutrino spectrum implies that 
that / will vary uniformly over the region [0,1]. 5 In that case the observed energy at the 
present era is 

E = I^- 

1 + 2 

fe 
l + z 

f -^unscattered (^*v 

where -E U nscattered = + z) would be the observed energy of the neutrino in the absence 
of resonant scattering. This shows how the entire allowed region of energy below E* is 
populated by rescattered neutrinos, with the energies shifted downward from the original 
spectrum. Especially interesting is the spectrum for i? U nscattered at or just below the limit 
E*: in that case, the only replenishment of flux is from the tail of the decay distribution 



We thank K. Hikasa for drawing our attention for this issue. 
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(/ ~ 1), from resonant production that has taken place only recently. (From Eq. (|T^, on can 
see that the conditions for this are that z ~ 0.) The restriction to / ~ 1 implies very little 
replenishment, so that a dip at E = E* should be a universal feature of the final spectra 
observed. This is a qualitative response, in the affirmative, to Question 2: the spectrum 
with absorption will show a dip at E — E*, and will be shifted downward from the spectrum 
absent resonant absorption. The complete effect of neutrinos emitted with non-resonant 
energies, passing through resonance, and then replenishing the flux at lower energies, is 
what we call accumulative resonance. 

4- A note on thermal broadening 

In the presence of the cosmological expansion, the effect of thermal broadening (because 
the CBN are not at rest) on our principal result (the universal dip described in the previous 
subsection) is negligible. The argument is as follows: after decoupling, the CBN spectrum 
is Fermi-Dirac, but in momentum rather than energy, even into the non-relativistic region. 
Thus the effect of thermal broadening is the introduction of a momentum spread in the 
target neutrinos of 0(T V ). The principal feature of our result, the dip at neutrino energy E*, 
occurs with neutrinos undergoing resonant scattering in the recent era, where the CBN is 
completely non-relativistic . In that case, the fractional energy shift of the target neutrinos 
is 0(T^/m1) ~ 10~ 3 (unless the lightest neutrino has mass < 10~ 4 eV). As can be seen 
from Eq. ([l3|) , the effect of this uncertainty is the introduction of a spread of (9(10 -3 ) in the 
value of z, the red shift at resonant scattering. The consequence, after integrating over red 
shift of the source, is that the sharpness of the dip at E* in the observed energy spectrum 
is softened by effects of O(l0~ 3 ) rather than 0(yl) ~ 10~ 15 , corresponding to the intrinsic 
width of the resonance. Thus, even with thermal broadening, the relative sharpness of the 
dip is preserved. 

5. Event Rates 

Next we consider the effect of the accumulative resonance on the total SRN differential 
flux. We first present the standard expressions for the SRN flux. Then we shall discuss how 
to incorporate the accumulative resonance effect. The differential flux of Supernova Relic 
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Neutrinos (SRN) is given by 

dF 
dE = 



ldN{t 



(1+*) 



dt 



dz 



dz , 



(20) 



(21) 



\ ^ 6 / e=(l+z)E 

where for heuristic purposes we adopt as our standard the Fermi-Dirac distribution 

dN(e) 120 e 2 1 

x x x exp (^y+T ' 

The constant £ = 0.5 x 10 53 ergs is the total energy carried by each flavor of neutrino. 
The temperature for the electron antineutrinos is Tp e N = 5 MeV and for the non-electron 
neutrinos and antineutrinos is = 8 MeV. However, a more general form 

'l + a) 1+a > 



dN 
de 



-(l+a)e/e 



(22) 



T(l + a)e 2 

has been proposed |24] which provides a good fit to simulated explosions [p5] of high-mass 
progenitors. Here where e is the average antineutrino energy at the source and the values of 



the fitting parameters e and a for the u e and u x spectra from three different groups [24, 25, 26 



(designated as LL, TBP, and KRJ, respectively), summarized in Table 1 of Ref. [27]. In 
Fig. [3] we will show the spread in our results obtained from the different spectra. 

Since we are considering the case of Dirac neutrinos, we take the CBN to consist of equal 
mixtures of left- and right-handed non-relativistic neutrinos and antineutrinos, with total 
number equal to two degrees of freedom in equilibrium during BBN. The resonant scattering 
will take place between right-handed SRN antineutrinos and CBN neutrinos, as well as left- 
handed SRN antineutrinos and CBN neutrinos. The problem becomes complex since the 
spectrum of SRN neutrinos and antineutrinos are different. Again, for illustrative purposes, 
we present results for the simplified case in which only SRN right-handed antineutrinos 
undergo resonant scattering, but from a CBN left-handed neutrino population given by n^, 
as defined after Eq. fll3f) . 

The Jacobian factor in Eq. )(|2~0D is given by 

km 1 1 



dt 

dz 



100- 



h(i + z)^n M {i + zf 



(23) 



s Mpc 

with Qm = 0.3, Qa = 0.7 and h = 0.7. The (comoving) rate of supernova formation i?sN is 
parameterized as 



^sn(^) 



0.013 



P* {z) 



(24) 
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where 

p, ( z ) = (1 - 2) x lO^Afoyr^Mpc" 3 x (1 + zf. (25) 

and the exponent changes at z = 1, from (3 ~ 2 — 4 (for 0<z<l)to/3~0 (for 
z > 1) The uncertainty in parameters describing Rsn(z) comes from the uncertainty in 
present knowledge of the Cosmic Star Formation Rate (CSFR) ||28|| . In this paper we choose 
"median" values for these parameters [p8| , Rsn(O) = 2 x 10 _4 yr _1 Mpc -3 , [3 = 2 (for 



< z < 1) and (3 = (for z > 1). 

The fact that the SN density is either constant (for far ones) or increasing with the 
distance (for near ones) is of great importance. It implies that most of the incoming neutrinos 
originate from distant SN (for far SN, the flux decreases like square of the distance while 
their density grows like the cube of the distance). Thus these neutrinos are redshifted and go 
through the resonance. The detection process is sensitive only to the flux of incoming anti- 
neutrino electrons. Consequently, to include contributions to the flux from the muon/tau and 
electron neutrinos we use the corresponding temperature in the expression for the neutrino 
spectrum, and include a mixing factor which for the electron neutrinos is 0.69 and for 
the muon neutrinos is 0.31. As discussed in Refs. [29, 3C, 3l]|, the relation between the v e 



spectrum observed on Earth to the various neutrino spectra at production depends critically 
on whether the neutrino mass hierarchy is normal or inverted. If normal (m 3 > m 2 > mi), 
then strong matter effects cause the v e at production to emerge from the stellar surface as 
the lightest eigenstate u\, with electron component \U e i\ 2 ~ 0.69. The small mixing of the 
electron with the third eigenstate \U e 3\ 2 <C 1 allows an equivalent two-flavor picture, with 
the result that neutrinos produced in the supernova as ^, x = fi or r will be received at 
Earth as v e with probability 0.31, and with energies corresponding to the v^jf spectrum at 
production. For the case of the inverted hierarchy, z/ e 's produced in the supernova emerge 
as the lightest mass eigenstate , now z/3. For sin 2 2#i3< 10~ 6 the resonance is non-adiabatic 
and there is complete conversion z/3 — > V\. This case then is the same as for the normal 
hierarchy. The adiabatic case (sin 2 2# 13 > 10~ 4 ) is very different: the original z/ e 's remain as 
P3 when emerging from the stellar surface, contributing negligibly to the v e flux at Earth. 
The entire v e flux at Earth then corresponds to the original v x produced in the supernova. 
For intermediate values of sin 2 2^13, the situation is of course more complicated. In this 
paper we will consider only the normal hierarchy. This can be regarded as conservative, 
since in some cases the v x spectrum (which generates the v e signal in the adiabatic inverted 
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hierarchy scenario) is harder than the v e spectrum, and will give rise to more events which 
escape the low energy cuts. 

In this context, we also note that Earth matter effects have been shown to modify the 
observed fluxes and spectra on earth ||32| . Since however the hierarchy in average energies 
between the v e and the other flavor is milder than thought this effect is expected to be 
subdominant Thus for simplicity and since these effects are not expected to induce 

gross modification of the observed spectrum we neglected them altogether. 

In order to obtain the observed spectrum, we note again that all neutrinos emitted from 
a source at redshift z with energies e outside the window E* < e < (1 + z)E* will arrive 
at z = without undergoing resonance, and with a flux given by Eq. (p0|) with energies 
E > E* and E < E*/(l + z). For source energies in the resonant absorption region, all of 
the original flux will undergo resonance absorption followed by decay into a flat spectrum. 
Since all energies, both before and after absorption, are redshifted at the same rate, one 
can obtain the rescattered spectrum by generating neutrino numbers according to Eq. fl2T|) 
in energy slices Ae at the source (for e in the absorption region), and redistributing this 
number according to Eq. (fi9| ) uniformly over the observed energy region < E < e/(l + z). 

In Fig. H] we show the resulting differential flux (with source flux given by Eq. (|21|) , with 
and without the accumulative resonant effect, integrated over redshift up to z = 4, and for 
Mq = 1.1 keV. As discussed following Eq. (|19D, there is a sharp dip at E = E* = MQ/2m v 
for all values of z. To demonstrate the spread introduced through differing assumptions 
about the source flux, we show in Fig. |] the flux with accumulative resonant effect for the 
three choices discussed in the context of Eq. (f22[). 

To estimate the event rates for SK and GADZOOKS |34| (SK enriched by Gd) we show 
in Fig. |] the differential neutrino flux folded with the detection cross section. This is 
for the inverse beta decay induced by the anti- neutrino capture in the detector |36[|.(For 
calculating the event rate we have used the quasielastic neutrino-nucleon cross section given 
in Ref. f37| .) The shape of the differential rate is modified due to the energy dependence of 
the cross section, which increases quadratically with energy. The main features of the effect 
due to the accumulation resonance such as the location of the dip and its width remain 
unmodified. 

The differential rate for SK and Kamland is rather low regardless of the presence of 
the accumulative resonance. We therefore present in Figs. the integrated flux for SK 
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and KamLand, with and without accumulative resonance, respectively. Since the energy 
thresholds for SK and KamLand are 18 and 6 MeV respectively, we note the interesting 
feature that the event rate for SK is roughly unmodified (the 18 MeV threshold is well 
above the absorption region) while a suppression of roughly 25% is obtained for KamLand. 
This feature gains definition in Fig. |j. The GADZOOKS experiment has a much lower 
background (due to the ability of identifying the emitted neutron) and therefore may be 
able to provide a differential rate information. 

dF/dE v (cm -2 s _1 MeV _1 ) 




FIG. 2: Depletion in the incoming SRN flux due to the resonance (solid curve) [j35|, compared to 
SRN flux without the resonance (dashed curve). The source flux is Fermi-Dirac. 

B. Non-resonance 

The effectiveness of the resonance process in redistributing the Rsn flux requires the 
resonance energy to be rather close to the peak energy E mSuX ~ 3Xf N of the Rsn (perhaps 
redshifted by a factor of 2) so that the boson mass must be in the range of 1 keV, 

A4_ = V*M1.57F) ~ 1 keV (^ ^ . (26) 
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We note that the BBN constraint ([!]) implies that the mass of the symmetry breaking 
scalar ~ lOkeV is typically above the resonance mass (not that far though). Thus it 
is more likely that the other light bosons (Goldstones or gauge) whose masses are almost 
unconstrained are required to provide a resonant channel for the Rsn scattering. 

In view of this restriction, it is important to check whether other non-resonant processes 
can become important. One interesting possibility within the present dynamical framework 
is shown in Fig. |7]. Again, this presumes the existence of either light Goldstone bosons \T2, 
IrJ|, |TT|, [TB| or light gauge bosons |14||. For s <C M? ~ f 2 we can estimate the cross section 



to be 

_ Uvf Vv Vv (r,<-r\ 

ml f 2 ml 

where we assume that the light bosons are produced on shell. This implies 



Mq < sj2m v Ej1 . (28) 
To have a substantial scattering we require the non-resonant mean free path A non _ res to be 




FIG. 3: Depleted incoming SRN flux due to the resonance for three different assumptions about the 
source flux. Dotted, solid and dashed curves are fits designated as LL, TBP and KRJ, respectively, 
after Eq. (||). 
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smaller than H 



-i 



Anon-rcs# ^ f < 1 , (29) 



which yields a lower bound on y v 



fc> l"^J ~ 10 tevJ ' (30) 

This requirement is valid for any value of Mq ■ 

The above process can have an effect on the SRN flux. If Mq < 2m u and there is sufficient 
optical depth, all the SRN will be transformed into invisible Goldstones and the signal is lost 
(for a related effect, see p8f). If Mq > 2m l/ then the process can effectively be characterized 
as v — > 4u, implying a substantial shifting of the entire SRN spectrum to lower energies. For 



a point source at distance l 1 the condition analogous to Eq. (30) for sufficient optical depth 



is 



fi /3000Mpc\ 1/4 . x 

2/,>3.3x 1CT 6 — v —\ (31) 




FIG. 4: Depletion in the incoming SRN flux folded with the cross section for detection with SK 
and GADZOOKS (dashed curve) compared to no resonance case (solid curve). Source flux is 
Fermi-Dirac. The essential features of the accumulation resonance remain unmodified. 
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FIG. 5: Integrated event rates for GADZOOKS and SK. The threshold for GADZOOKS is 10 MeV 
and for SK it is 18 MeV. Source flux is Fermi-Dirac. 
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FIG. 6: Integrated event rates for KamLand. The threshold for KamLand is 6 MeV. 
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FIG. 7: Two Goldstones production through Higgs exchange, off-resonance process. 



where I is the distance travelled by the SRN. For SN1987A, £ = 50, 000 pc, the fact that non- 
resonant scattering have not occurred, i.e. neutrinos with undegraded energy were observed 
3~9fl , gives an independent upper bound on y u , 



5.5 x 10~ 5 . (32) 

This bound is comparable to the cosmological one in Eq. QT]) . However, considerably more 
detailed work is required in order to establish such a bound: a combined likelihood analysis in 
the symmetry breaking scale and the parameters describing the neutrino spectrum needs to 
be done in order to establish confidence levels for all variables [ 40 1 , followed by marginalizing 



on the spectrum parameters. In the meantime, we adopt (|32"D as a rough, and preliminary, 
indication that this bound can be comparable to others we have mentioned. 



III. CONCLUSIONS 



In the models in which neutrino masses are light due to the low-energy symmetry breaking 
scale, extra light bosons are typically present. These light bosons couple to neutrinos 
with the coupling that is proportional to their masses and therefore directly related to the 
symmetry breaking scale. We have shown that, in principle, one can measure this coupling 
because SN neutrinos interact with cosmic background neutrinos via these bosons modifying 
the SN neutrino flux dramatically. 

We have discussed two types of processes that are present due to these interactions. The 
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first is a low energy analog of a Z-burst [[ITJ where neutrinos interact producing an on-shell 
boson which subsequently decays to a pair of neutrinos. The expansion of the universe allows 
for a wide range of energies in which such a mini Z burst can occur. We characterize this 
process as accumulative resonance. The second is a non-resonance process which leads to a 
global degradation of energies in the supernova neutrinos flux. The observation of neutrinos 
from the 1987a supernova yields an important constraint on the parameter space of the 
above models since the observations were fully consistent with no such degradation. 

As is often the case with observations related to neutrino physics, the signal we find is 
currently beyond the limit of each individual present experiment. The strength of our signal 
can improve once data from several experiments are combined. For example, our analysis 
reveals that a robust prediction of these models is that results from SK will be unaffected 
by the above processes while KamLand should observe a suppressed flux. However, in order 
to be convinced that depletion is observed it is desirable to actually observe the predicted 
dip in the flux, which requires certain amount of information on the energy dependence of 
the flux. This can be obtained in the future via water Cerenkov detectors enriched with 
added Gadolinium, (the GADZOOKS proposal). This can be implemented in the very 
near future using an upgrade of SK, or attained via future experiments such as HyperK, 
UNO etc... Such and experiments can collect tens of relic supernova neutrinos per year and 
provide us with information on the energy dependence of SN flux. Furthermore they are 
more sensitive to observation of a single SN event. In principle these Megaton scale water 
Cherenkov detectors might detect neutrino burst from O (Mpc) distant SN with average 



rate of one burst per two years or so (see e.g. [|2| and references therein). In this case there 
are two possible scenarios. The first, less exciting, is that the future burst of SN neutrinos 
will have differential flux consistent with the one observed from SN1987A. In this case we 
expect the bound on the model parameters to be improved. Alternatively, a more exciting 
possibility is that the incoming neutrinos are absent, or are severely degraded in energy 
relative to expectations. This can then be interpreted as a measurement of the coupling 
between the neutrino, the Higgs and the Goldstone bosons which yielded the degradation 
in the flux through the non-resonance process. 

We note parenthetically that we have for simplicity neglected modification of the neutrino 



spectra due to shock wave effects fI5[] . These effects can change the spectral features, 
inducing non-adiabatic transitions, in a manner that depends on the neutrino flavor 
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parameters [44|. However, since the exact dynamics related to the shock propagation in 
the SN is not well understood and since the above effects are red-shift dependent and will 
be smeared out once the integration over z is applied we shall not include then in our 
computations. In this context it is important to emphasize that the position of the universal 
dip will not vary with redshift. 

It is also important to stress that in this work we focus on the modification of the dynamics 
of the SN neutrinos while they propagate in space outside the SN. Interesting effects may 
be induced by the presence of the new light degrees of freedom inside the SN [55] . These, 
however, are model dependent; for example, in the case of late Dirac neutrino masses the 
overall effect is expected to be miniscule |12|]. The Majorana case is more involved since 
the extra bosons are expected to be thermalized inside the SN core fl2| , |45| through their 
reactions vv <-> G. Consequently, the dynamics inside the SN might be significantly modified 
due to the presence of new flavor and lepton violating interactions. This requires a more 
detailed study which is beyond the scope of our paper. Note that it is likely that only 
the (pseudo) Goldstone boson will be light enough to be thermally produced inside the 
core [|1^, [L3| . Thus the modification of the the neutrino spectra [|46| is expected to be below 
the 10% level expected just by counting degrees of freedom, which is probably smaller than 
other systematic effects which have not been included in this analysis. 

A summary of the available parameter space, subject to BBN bounds, for observation of 
the resonant process is given in Fig. [| It is clear that BBN constraints are most constrictive 
in the case of Majorana neutrinos with a non-Abelian symmetry breaking sector. Least 
restrictive is the Dirac, Abelian, scenario; it also has the least impact on SN dynamics at 
the source. 

We stress that observation of the above signals may not only shed light on the origin 
of neutrino masses but also yield an indirect observation of the elusive background relic 
neutrinos. In addition, we note that the above processes may be induced by other light 
degrees of freedom which couples to neutrinos. Consequently our signal may provide a test 
for other frameworks apart from the late neutrino masses one. 
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FIG. 8: The cosmological bounds and the regions for the supernova neutrino spectrum distortion 
due to the resonance and non-resonance processes for a single Major ana (Dirac) neutrino for 
an abelian (non abelian) model are shown in {jj u ,Mq) plane. The region above the red (solid 
horizontal) line is excluded by the BBN constraint (for the Dirac case), SN cooling (for Majorana 
case) and due to the observation of (undegraded) SN1987A neutrinos. In the region below the blue 
(solid slanting) line the mean free path is too long for the resonance to occur. The region above 
the green (dashed slanting) line, which is relevant only for the non abelian Majorana case, is the 
region excluded by the BBN constraint. The region above the black (dashed horizontal) line is 
the region of the future experimental sensitivity to the observation/non-observation of the SRN 
neutrinos due to the the non-resonant processes. 
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